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ABSTRACT 


I 

'~'We  examine  a  quantum  mechanical  model  for  the  enhancement 
of  fluorescence  yield  and  photochemical  rates  by  the  presence  of 
small  particle  capable  to  sustain  electromagentic  resonances. 

The  sphere  lowers  the  absorption  and  emission  efficiency  by 
taking  energy  from  the  molecule  and  storing  it  into  nonradiative 
modes,  enhances  the  absorption  by  increasing  the  local  field  and 
enhances  the  fluorescence  by  emitting  efficiently  the  energy 
transferred  from  the  emitting  state  to  a  radiative 
electromagnetic  resonance  of  the  particle.  The  present  paper 
studies  how  the  Interplay  of  these  effects  modifies  fluorescence 
intensity  and  the  rate  of  photo-chemical  reactions^ 


1. 


INTRODUCTION 


The  spectroscopy  and  the  relaxation  dynamics  of  vibrational 

and  electronic  excitations  of  molecules  and  atoms  adsorbed  on 

metal  surfaces  have  received  a  great  deal  of  attention  lately 

since  they  are  involved  in  many  surface  processes  such  as  surface 
1“4  5 

photochemistry  ,  surface  photoemission  ,  surface  enhanced  Raman 
scattering  (SERS),®'^  and  photolumin^scence . In  the  present 
article  we  examine  molecular  fluorescence  and  photochemistry  in 
the  presence  of  a  surface. 

Qualitatively  these  phenomena  are  fairly  well  understood. 

In  all  surface  enhanced  processes  the  local  field  acting  on  the 

molecule  is  increased  by  reflection  of  the  incident  radiation. 

This  increase  is  substantial  if  surface  curvature  is  large, if 

the  incident  light  excites  surface  electromagnetic 
1 1  —  1. 6 

resonances,  or  if  solid  particles  are  used  to  polarize  each 

17 

other.  This  local  field  enhancement  increases  the  rate  of 
photochemical  processes  and  the  intensity  of  adsorbate 
fluorescence.  The  presence  of  the  surface  however  causes  energy 

transfer  from  the  molecule  to  surface  excitations,  such  as 

18^13  20—21 

plasraons,  electron-hole  pairs  or  phonons.  This 

results  in  a  decrease  of  the  excited  population  which  causes  a 

decrease  of  the  photo-chemical  rate.  The  total  rate  is 

established  through  the  competition  of  these  opposing  effects. 

In  the  case  of  fluorescence  an  additional  complication 

arises  because  the  excited  molecule  polarizes  the 
7a  8  3  18  13 

surface  •  •  >  >  and  this  emits  at  the  molecular  frequency, 

thus  enhancing  the  emission  intensity.  This  enhancement  is 
particularly  strong  if  the  emission  frequency  matches  that  of  a 
radiative  surface  resonance. 

Recently  several  theoretical  models^  '  predicting 

Q  — 1 0  •  24 

and/or  explaining  observed  fluorescence  'and  photochemistry 


of  adsorbates  have  appeared  in  the  literature. 


Our  purpose  in  this  article  is  twofold:  to  provide  a  more 
systematic  derivation  of  the  quantum  mechanical  behavior  of  the 
molecule-sphere  system,  and  to  carry  out  a  more  extensive 
numerical  study  of  fluorescence  and  photochemical  rates. 

Previous  work  included  surface  effects  in  an  intuitive  manner  by 
replacing  the  oscillator  strength  and  width  appearing  in  the 
quantum  theory  for  an  isolated  molecule  with  surface  dependent 
expressions  suggested  by  the  classical  theory.  The  intuitive 
arguments  are  quite  reasonable,  but  we  prefer  to  derive  the  rate 
equations  by  using  a  detailed  quantum  model.  Our  results  confirm 
those  of  previous  work. 

The  theory  uses  a  small  (Rayleigh  limit)  sphere  as  a  model 

of  the  surface.  This  can  be  experimentally  realized  in  colloidal 

2  S  2  6“ 2  7 

systems  or  by  preparing  matrix  isolated  metal  spheres.  To 

some  extent  a  small  sphere  is  a  crude  model  of  a  rough  surface; 

however  such  extrapolations  must  be  made  with  care  since  all  the 

effects  mentioned  above  could  depend  strongly  on  the  shape  of  the 

surface. 


11. 


THE  MODEL 


We  develop  a  quantum  mechanical  model  for  the  steady  state 
fluorescence  or  pho.tochemistry  of  molecules  adsorbed  on  a  metal 
substrate.  The  system  we  have  in  mind  consists  of  a  small 
metallic  particle  or  a  rough  metal  surface  covered  with 
adsorbates  and  irradiated  with  an  external  light  source.  We  are 
concerned  with  the  influence  of  the  substrate  on  photochemistry 
of  adsorbates.  To  this  end  we  consider  the  incoherent  light 
scattering  from  a  molecule,  idealiied  as  a  three  level  system, 
adsorbed  on  a  dielectric  substrate. 

II . 1  The  Three-Level  System 

The  schematic  diagram  of  the  adsorbate-substrate  system  is 
shown  in  Fig.  1.  State  | I>  is  the  initial  state  (ground 
electronic  level)  and  the  states  |II>  and  |II1>  are  excited 
electronic  states.  These  are  molecular  states  "dressed"  by  the 
presence  of  the  substrate. 

State  I II>  is  populated  by  absorption  of  energy  from  the 
light  source.  This  is  accompanied  by  a  radiationless  relaxation 
to  state  |III>  and  radiative  relaxations  to  state  |I>.  If  state 
|III>  is  a  dissociative  continuum  the  molecule  undergoes 
photochemical  decomposition.  If  1III>  is  a  bound,  discrete 
electronic  level  it  can  decay  radiatlvely  to  the  ground  state 
|I>.  The  presence  of  the  surface  modifies  the  transition  |I>  -*■ 
|II>  and  |II>  1III>  but  leaves  11I>  -*  |I1I>  unchanged.  The 
latter  assximption  excludes  the  chemical  effects  of  the  surface 
and  it  might  be  inaccurate  for  molecules  chemisorbed  at  the 
surface.  This  is  not  a  serious  limitation  here  since  we  are 
concerned  with  long  range  electromagnetic  effects  taking  place 
many  A  away  from  the  surface . 


II. 2  Absorbing  System:  Spontaneotis  Emission 


We  consider  the  Interaction  of  the  external  radiation  with 
the  molecule  (only  absorbing  system  |I>  |II>  is  considered)  to 
determine  the  rate  of  absorption  and  the  induced  molecular 
dipole.  The  latter  must  be  determined  by  a  self  consistent 
procedure  which  includes  self-polarization  caused  by  the  fact 
that  the  induced  dipole  polarizes  the  surface  and  interacts  with 
that  polarization  field. 

The  incident  field  ^^(u,t)  is  reflected  by  the  surface  and 

this  produces  a  new  field  denoted  here  u , t ) ( t , w ) ,  and 

called  the  reflection  field.  A  large  number  of  formulae  for  the 

reflection  tensor  "^(u)  for  various  surface  geometries  are  given 

in  Reference  7a.  The  molecular  dipole  induced  by  ^^(u,t) 

polarizes  the  surface  and  the  polarization  field  at  the  point 

7a  1 4 

where  the  molecule  is  located  is  denoted  ' 


^  (w,t)  •  <li(t,«)> 

"^(u.t)  is  called  here  the  image  tensor  and 

<l!(t,u)>  2  \  []I_(w)e"^“^  +  C.C.]  (II. 1) 


is  the  Expectation  value  of  the  induced  molecular  dipole.  To 
calculate  the  Induced  dipole  self-consistently  we  include  the 
field  u ,  t )  •  <li(  t ,  u )  >  in  the  expression  of  the  local  field 


+  T*  (u)  •  <M  (u)>]e  c.c. 


(II.  2) 


which  is  responsible  for  the  polarization  of  the  molecule.  In 
this  expression  I  is  the- unit  tensor.  The  Hamiltonian  therefore 
contains  a  term  of  the  form 


2TTh«^  -=  . 

Vj(t)  =  {i  Z  (-  y  [12><llM2i-(t'+  R'((Ojj))-e-  '“21’^ 

ic 

-  Ic  12><ll|l2i-^(u)*<^_{u)  >e“^^“'“21^^+  c.c.]}  .  (11.3) 


where  V  is  the  volume  of  the  .pavlty  used  for  field  quantization; 

the  svim  is  carried  over  all  the  photon  modes  in  the  cavity;  a^  is 

the  annihilation  operator  for  a  photon  of  frequency  and 

polarization  denoted  by  a;  is  the  unit  polarization  vector. 

is  the  transition  dipole  moment  corresponding  to  a  molecular 

transition  |i>  -*>  |j>,  and  =  (Ej-E^)/h.  In  writing  Vj(t)  (see 

Equation  II. 3),  which  is  the  molecule  field  coupling  in  the 

interaction  representation,  we  use  the  rotating  wave 

approximation  which  neglects  small,  noh-linear  terms.  The  image 

field  -has  been  treated  semi-classically  while  the  incident  and 

the  reflected  fields  are  quantized.  The  quantization  is  needed 

only  to  make  sure  that  spontaneous  emission  is  treated  properly; 

other  quantum  electrodynaraic  effects,  are  not  expected  to  play  a 

role  in  surface  spectroscopy  where  working  at  intensities  where 

the  statistical  properties  of  the  photon  become  important  is  not 

practical.  The  "image"  term  was  treated  semiclassically  for 

convenience,  since  otherwise  we  would  have  to  quantize  the 

electrons  and  phonons  in  the  solid  to  get  the  correct  expression 

for  the  mariner  in  which  the  polarized  solid  renormalizes  the 
2 1 

photon  field.  Note  that  in  treating  the  semlclassital  term  we 
must  disregard  the  photon  states  in  the  wave  function.  Thus  in 
dealing  with  the  quantized  field  we  have  ]!>  =»  Il>|nj^+1>  and  |II> 
a  |2>|nj^>  where  ]  1>  and  12>  are  the  molecular  wave  functions  and 
|nj^+l>  and  lnj^>  are  field  wave  functions  in  the  occupation  number 
representation  (a  coherent  state  would  be  more  adequate  in 
describing  laser  radiation).  However,  in  dealing  with  the 


semiclassical  term  ] I>  =  |1>  and  |II>  =  |2>. 


The  equations  of  motion  for  the  matrix  elements  of  the 


density  operator  are 


Pii{t)  =  -  J  <il[V(t).p(t)l|i>  .  I  , 


(II.4) 


p^^(t)  =  (-iu.j-r.j)p,j(t)  -  J  <iitv(t)  ,p{t)]  1  j>,  i  pi  j 

(II. 5) 

where  {ii>,  lj>}  =  {II>.  |II>>.  V(t)  is  the  'interaction 
potential  given  by  (II. 3),  but  taken  in  Schrodinger 
representation;  is  a  phenomenological  parameter,  that 

characterizes  depopulation  induced  by  interaction  with  the 
medium.  In  the  interaction  representation,  Eqs .  { II . 4 ) - ( II . 5 ) 
can  be  written  as 


pJjCt)  = 


(II. 6a) 


(II. 6b) 


Since  <p(t)>  =  Tr(pp(t)),  the  induced  dipole  <|i(t)>  depends  on 
itself  through  p(t).  However,  Eqs.  (11.5),  (11.6)  and  (11.3)  can 
b'^  solved  to  obtain  self-consistent  expressions  for 
P2^(u)  and  <M^{w)>: 


,1  .  U_.»-  G  (.1  . 

^  -P.j)  ; 


(II. 7a) 
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(II. 7b] 


<P_(u)>  (Pll"p22^  >in+l  [^21-  (^+^(1^)  )-e] 


[l(tj-w^^+ir^^)  + 


^^11~^22^ 


(II. 8a) 


^12"21-'®<“>'''''^12 


<p^((j)>  =  <l5_(u)> 


(11.8b) 


In  these  expressions  the  intensity  of  the  light  source 
appears  through  the  factor  n+i,  where  n  is  the  photon  occupation 
number  of  the  incident  field,  which  is  given  by 


n+1  =  V|E^|  ^(STThu)"^ 


(11.8c) 


The  term.R(u)  is  particularly  large  if  the  light  excites  a 

radiative  electromagnetic  resonance  of  the  surface.  The 

polarization  of  the  surface  by  the  induced  molecular  dipole 

appears  through 'g(u) .  This  quantity  acts  as  a  molecular  self- 
21  3  Os 

energy  '  since,  by  polarizing  the  surface  and  interacting  with 
the  polarization  charge,  the  molecule  polarizes  itself.  As  in 
any  other  areas  of  many  body  theory  the  self-energy  shifts  and 
broadens  the  original  state.  Note  that  even  though  we  treated 
the  sphere  by  classical  electrodynamics  the  structure  of  the 
theory  is  the  same  as  the  one  obtained  by  a  full  quantum 
treatment  of  the  electrons  in  the  sphere'  and  the  two  level 
system.  The  numerical  values  of  the  self -energy  are  however 
different:  a  full  quantum  theory  for  the  sphere  introduces  non¬ 
local  electromagnetic  effects,  which  are  important  at  small  (i.e. 
30 

about  5A  )  distances.  A  detailed  comparison  between  local  and 
non-local  theories  can  be  found  in  Ref.  30b.  Throughout  this 
paper  we  use  for  Giu)  a  local  theory.  In  this  theory  "g(«)  is 


enhanced  whenever  the  dipole  couples  to  the  radiative  or  the  non- 
radiative  resonances  of  the  surface.  Note  that  if"^  = *R  =  0  the 
expressions  derived  above  reduce  to  those  for  a  two  level  system 
in  the  absence  of  the  sphere.  The  width  T contains  the 
dephasing  and/or  the  depopulation  caused  by  intramolecular  and/or 
radiative  transitions  and  the  interaction  with  the  medium. 


It  is  interesting  to  note' that  the  term 


{(u-U2i  +  ir2i)  (Pii-p22^/^^^12^2r 
is  a  3  X  3  tensor  since  the  symbol 

^12^21’  is  a  3  X  3  tensor  whose  xy  component  is 


^^21^x°xy  ^^21^y°xy‘^^‘'‘21^z°zy^ 


Thi'S  is  introduced  in  the  theory  by  the  presence  of  the  image 
field  in  the  self-consistent  calculation  of  <p>  and  P2i' 


tensofial  character  of  this  quantity  does  not  mean  that  the  life¬ 
time  is  tensor.  If  one  performs  all  the  tensor  Inversion  and 
multiplication  operations  required  in  (II. 7a)  the  consequence  of 
this  formula  is  that  the  life-time  depends  on  the  orientation  of 
the  transition  dipole  with  respect  to  the  surface.  In 
particular,  parallel  and  perpendicular  dipoles  have  different 
life-times.  Since  the  life-time  depends  on  the  rate  of  energy 
transfer  from  the  dipole  to  the  surface,  the  existence  of  the 
above  effect  in  the  theory  is  reasonable. 


The  induced  dipole  <u_(u)>  and  the  off-diagonal  element  P2, 
depend  on  the  population  difference  (Pii~P22^*  When  the' 
populations  of  the  ground  and  excited  states  become  equal,  the 
induced  dipole  moment  of  the  molecule  vanishes  and  so  does  the 
surface  induced  broadening  and  frequency  shift.  This  is  the  well 
known  saturation  effect  caused  by  intense  radiation.  If  the 


external  field  is  weak  enough,  W:e  can  approximate  - 

o  .o  ,  _  o  ^  o 


~  ^22^  where  s:  1  and  n  0  are  thermal  eouilibrium 


values  of  the  populations.  We  shall  not  resort  to  this 
approximation  since  the  strength  of  the  surface  enhanced  local 
field  can  be  large  for  certain  choices  of  parameters. 


The  above  calculations  provide  us  with  self-consistent 
expressions  for  <n>  and  compute  the  rate  of  change  of 


population  we.  take  the  external  field  to  be  [T  +  ^  (u)]"?^(u) 


This  expression  does  not  contain  the  field  G"<p>  since  the  latter 
is  already  included,  by  the  use  of  the  self-consistent  procedure, 
in  the  self-energy  of  the  two  level  system.  Using  the  expression 
of  Vj{t),  modified  as  indicated  above,  and  Eqs.  (11.4)  and  (II. 7) 
we  obtain 


CPjs^^Jlrad  "  (n+1)  (4TT«/hV)  (p^^-P22)  iPgl 

[(l-?i“^li2j-‘^(«)  (II.  9) 


and  Pj^^  (t)  =  ”P22  gives  the  the  rate  of  change 


of  the  population  of  state  (II>  by  the  emission  of  a  photon  of 
frequency  u,  polarization  e  and  wave-vector  In  addition  to 


radiative  contributions  to  P22^^^  there  will  be  non-radiative 


decay  channels  which  are  not  Included  in  Eq.  (II. 9). 


Although  not  necessary  in  principle,  we  shall  make  certain 
simplifications  in  the  mathematics  assuming  that  both  the  induced 
dipole  and  the  incident  field  have  the  form  p  »  p  z  and 

•  m 

S  »  E  z,  where  z  is  the  unit  vector  along  the  z  axis.  Then  the 
o  o 

image-field  due  to  the  dipole  induced  in  the  surface,  at  the 
location  of  the  dipole,  picks  up  only  the  z-component  of  the 
induced  dipole.  So  only  G  (u)  s  G(u>  is  non-zero. 


10 


Thus  we  can  rewrite  the  equations  derived  here  .as 


P2i(«)  =1  [ii2i- {T+^(u)  )-fe] 

<■?.{«)>  =  <PirP22^ 

(^21*  J Pl2'^^““"2l'*‘^’’21^ 


(11.10) 


(11.11) 


and 


■  '■safll  -  ''■se  '■sp’>'>22  ' 


rad 


where 


(11.12) 


■>'21^“^  *  ^21  ■  K  (Pii~P22’ ilniG(u)]  , 


(11.13) 


r  and  r  are  the  stimulated  absorption  and  emission  rates  for 
sa  se  j-j 

the  transition  |I>  •*-*  |1I>/  and  r"  is  the  corresponding 

sp 

spontansous  photon  emission  rate.  The  explicit  expressions  for 
various  rates  of  absorption  and  emission  are  obtained  by 
inspection  of  Eq.  (II. 9): 


sa 


r__  =  r®p(n+i)  , 

n  , 
sp 


se 


rl2>  ,iis, 


where 


y2i(«)  s  121^“^  *  ^'^21^“^ 


(11.14) 

(11.15) 

(11.16) 

(11.17) 


and  y2i(“)  121^“^  real.  For  large  n,  Pgg^-^'sg 

expected.  They  depend  on  the  incident  laser  intensity  through 
the  relation  n  “  ]  ^/8TTh« . 


The  total  spontaneous  emission  rate  is  obtained  by  summing 
Eq.  (11.16)  over  all  wave  vectors  jc  and  polarizations  a.  Using 
the  fact  that  the  number  of  modes  per  unit  frequency  interval  and 


per  xinit  solid  angle  is  (for  a  given  polarization) 
2 


d^No 


Vu' 


dwdn 


(2Tr)^c^ 


(11.19) 


we  obtain  the  spontaneous  emission  rate  per  unit  solid  angle 
(about  the  propagation  direction)  and  per  unit  frequency  range: 


d^r^^^  3 

sp  _  2(i» 


lli2i*(r+r(«))*e|2 


dudO 


Tihc^  t(<*J-U2j^-y22^(«)  )^+r2i(“)^] 


(11.20) 


In  order  to  get  the  total  spontaneous  emission  rate  "^we 


sp 


integrate  the  right  hand  side  of  Eq.  (11.20)  with  respect  to  u 
and  n.  Approximating,  as  customary,  1^21  ^“^“^21  ^“21  ^  '  where  U21 
is  the  transition  frequency,  the  total  spontaneous  emission  rate 
is 


j,(2)T^  _4_ 


sp 


he' 


‘'.3 

“21 


i;2i-{nr(C'2i))-  ^ 


(11.21) 


where  «2i  =  "21  ""  ’*21 


(W21);  721  shift  of  the  upper 


level  due  to  the  presence  of  the  surface.  This  shift  is  usually 

II 

small  compared  to  «2i  hence  may  be  neglected.  When  r2i 
neglected,  the  effect  of  the  sphere  on  the  spontaneous  emission 


(within  the  approximations  made  here)  appears  through  ^(“21^ 


we  recover  the  usual  formula  when  the  sphere  is  absent : 


.(2)T 


sp 


—  I  u  •  e  !  ^ 

3  “21 ' ^21  ® ' 
nc 


(11.22) 


which  is  identical  with  Einstein's  A-coef f icient .  A  similar 

inspection  of  Eq.  (11.15)  in  the  absence  of  the  surface  and  in 

3 1 

the  limit  r2j^-*0  gives  Einstein's  B-coeff icient  ,  implying  that 
the  presence  of  a  polatizable  object  renormalizes  the  Einstein 
coefficients  via  3^^  and  ( ^ ) ) . 


11.3  STEADY  STATE  POPULATION 


So  far  we  have  focused  our  attentioR  on  the  absorbing 
system  (|I>  •*-*•  |II>)  alone  and  established  expressions  for  the 
power  absorption  by  the  molecule  in  the  presence  of  the  sphere  as 
well  as  the  rate  of  spontaneous  emission  of  photons  via  the  |II> 

-►  |I>  transition.  These  formulae  contain  the  level  populations, 
which  we  now  set  out  to  determine.  The  emitting  system  consists 
of  levels  |III>  and.  |I>  and  fluorescence  emission  takes  place  via 
the  spontaneous  transition  |III>  -*■  1I>.  Since  the  fluorescence 
intensity  is  low  we  do  not  include  stimulated  processes  in  the 
|III>  |I>  transition.  The  absorbing  and  emitting  systems  are 
coupled  only  through  the  radiationless  transition  |II> 

We  denote  by  and  the  rates  of  decay  of  populations  of 

levels  |2>  and  |3>,  respectively,  for  an  isolated  molecule 

modeled  as  a  three  level  system.  They  Include  both  radiative  and 

nonradlative  parts,  such  that  we  can  write  r,  ■  Ti  ' 
i NR ]  ■  11 

•  (i  *  2,3),  and  treat  the  full  as  an  empirical  molecular 

parameter  of  the  theory.  The  presence  of  the  sphere  simply 

renormalizes  these  decay  rates,  by  adding  to  them  the  surface 

induced  terras  discussed  in  Section  II.  1.  We  therefore  denote 

their  values  when  the  sphere  is  present  by  and  r^.  Written 

out  explicitly  these  quantities  are 

^2  “  ^2  —  ^Pll"P22^  ImG(u)  {II. 31a) 


r3  -  r3  H.  - ^  (Pn-P33^  (ii.aib) 

where  u  is  the  frequency  of  the  incident  laser  and  u’  that  of  the 
fluorescent  emission.  Further,  let  us  write  r  =  r*  .  Then  the 

)  S  C  So 

equations  of  motion  for  the  populations  are 


(II. 32a) 
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^se^P22“Pll^ 

^2^22 

^3P33 

(II. 32a) 

P22 

“  ^se^Pll“P22* 

^2^22 

-  ^P22 

(II. 32b) 

P33 

KP22  -  f'sPaa 

(11.32c) 

steady  state  solutions  for  p.. 

(1  =  1.2,3) 

are  obtained  by 

requiring  that  =  P22 

“  P33 

=  0: 

Pll 

*  (rgg+K+r2)/[K+r 

2^^se 

(2+K/f3] 

(II. 33a) 

P22 

=  r  p, ,/(r  +K+r 

se^ll  '  se 

2^ 

(II. 33b) 

P33 

(  )  P22 

*  in 

( II . 33c) 

The  level  populations  are  dependent  on  the  molecule-surface 
separation  through  and  (i  =  2,3).  These  steady  state 
values  of  the  level  popu-ations  can  be  used  to  compute 
(Eq.  (11.13)). 


II. 4  ENHANCEMENT  OF  FLUORESCENCE 


The  spontaneous  radiative  decay  of  state  | III>  causes  a 
fluorescence  signal  of  frequency  The  polarization  dipole 

induced  in  the  sphere  by  the  emitting  molecular  dipole  radiates 
at  the  source  frequency  u'.  Thus  the  radiation  from  the  system 
consists  of  the  coherent  siim  of  the  radiation  for  both  molecular 
dipole  and  the  dipole  induced  in  the  sphere.  However,  instead  of 
going ’ through  the  complications  of  determining  the  molecular 
(emitting)  dipole  moment  at  u',  and  computing  the  fluorescence 
from  the  classical  radiation  formula,  we  write  the  power  emitted 
by  the  molecule  due  to  a  spontaneous  transition  |III>  -*•  1 1>  as 

P  =  Np  hu'  ,  (11.34) 

FI  33  sp 

where  N  is  the  total  number  of  molecules  in  the  system  and  f 

sp 

is  the  total  rate  of  decay  due  to  spontaneous  emission  from  state 
|II1>,  and  is  given  by 

(i:.35a) 

or 

(II. 35  b) 


In  writing  Eq.  (II. 35b)  we  neglected  the  frequency  shift  caused 
by  the  sphere  and  averaged  over  the  random  orientations  of  the 
polarization  direction  e.  If  we  denote  the  power  emitted  by 
isolated  molecules  by  P-,,  then  the  enhancement  ratio  for 
fluorescence  is 


FI 


FI 


*F1 


11"^I“31>|"<P33/P3?' 


(i:.36a) 


where 


(o) 

P33 


is  the  population  of  level  3  in  the  absence  of  the 


(o)  (o)  — 

sphere.  Since  (P33/P33)  *  (P22/p22^  can  write 
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2  2^3  ^21<“J  [(«-«2l)^^(<r'  )2] 

|1+R(«)|2  - j 

Tg  21  [(u-Wjj^)  +r2j(u)  ] 


K+r2+r^°^2+K/r3) 

K+r2+r^e^2'^^''^3^ 
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(II. 36b) 


In  a  r-esonance  situation  (u2«d  and  w'stw.,,)  we  neglect  the  level 

fR)  ^ 

shifts  and  assvune  r  <<  r'  '  ,  to  obtain 
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R_.S  |1+R((J)  |^|1+R(U')  (;;;^)  (— - 
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Here  «  (r2+K)/2 


( 11.36c) 
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Let  us  now  analyse  Eqs.  (11.36).  |ltR(u)|  is  strongly 

enhanced  when  the  incident  laser  field  is  in  resonance  with  a 

radiative  (i.e.  dipolar)  electromagnetic  resonance  of  the  sphere. 

This  implies  enhanced  absorption  and  consequently  enhanced 

emission.  Another  source  of  enhancement  of  emission  is  through 
2 

|l+R(w')|  ,  which  is  large  when  the  emission  frequency  is 
resonating  with  a  radiative  electromagnetic  resonance  of  the 
sphere.  If  the  separation  between  the  excitation  and  the 
fluorescence  frequency  is  much  smaller  than  the  width  of  the 
electromagnetic  resonance  a  double  resonance  condition  can  be 
achieved.  On  surfaces  having  more  than  one  electromagnetic 
resonance  (e.g.  ellipsoidal  particles)  one  can  try  to  use  one 
resonance  for  enhancing  excitation  efficiency  and  another  to 
enhance  emission. 


As  noted  earlier,  the  enhancement  is  limited  by  the 
broadening  of  the  upper  level  due  to  the  presence  of  the 


substrate.  Due  to  the  bradening  of  state  | II>  absorption  is 

I 

decreased  as  reflected  by  the  factor  j/ygj (“) ) •  Similarly,  the 
broadening  of  State  |  III>  brings  in  a  factor  (Tg/Tg)  which 
counteracts  in  part  the  enhancement  caused  by  |l+R(u')l  •  The 

remaining  factors  are  presumably  due  to  the  decrease  in  the 
transfer  of  population  from  ill>  to  |1II>.  Such  qualitative 
conclusions  can  be,  of  course,  anticipated  on  the  basis  of  simple 
classical  models. 


II. 5  Photochemical  reactions  and  radiationless  transitions  for 


adsorbates 


Photochemical  reactions  and  radiationless  transitions  are 
very  similar  processes  in  which  the  molecule  excited  in  state 
|II>  undergoes  a  transition  to  a  continuum  or  another  bound 
state,  respectively.  We  consider  here  cases  in  which  the 
presence  of  the  surface  does  not  affect  the  transition  rate,  thus 
excluding  catalytic  processes.  The  surface  affects  the  observed 
yield  by  modifying  the  primary  process,  i.e.  the  population  of 
the  level  |II>. 

It  is  therefore  sufficient  to  consider  a  two  level  model 

29 

with  the  equation  of  motion 

"ll")  =  ♦  "‘ae  *  ^  *  ''pc’'’22  ' 


”22"’  “  ■‘’11"'  •  (i:.37b) 

We  have  asstimed  that  the  reagent  concentration  is  high  enough  for 
the  photochemical  reaction  to  have  any  effect  on  the  steady  state 
kinetics,  so  that  -  1.  The  steady  state  solution  of 

Eq.  (11.37)  yields 
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{I-P22)  ' 


( II . 3aa) 

(II. 38b) 


where  f  ^  r  .  Since  K  is  the  number  of  times  a  molecule 
sa  se  pc 

undergoes  photochemical  reaction  per  second,  the  number  of 
molecules  decaying  per  second  by  photochemical  decomposition  is 
given  to  be  Ii'  units  of  N,  the  photochemioal  rate  is 


«  o 


(11.39) 


Vse 


The  enhancement  factor  for  this  rate  Is  the  ratio  R  =  r  / 

pc  pc 


Here  r2j^.  which  enters  through  and  Tj.  is 


Ill .  NUMERICAL  RESULTS 


We  shall  now  present  numerical  results  for  a  molecule 

located  at  a  distance  d  from  a  spherical  Ag  particle  of  radius  a. 

The  dielectric  constant  of  the  sphere  is  taken  from  Ref.  32.  The 

molecular  parameters  are  the  widths  r.  and  T-,  the  oscillator 

^  2  2  ^ 

strengths  f^^  and  f^^^  ^^ij  ~  ®  ™ 

electron  charge  and  mass  respectively) ,  and  the  rate  K  for  the 
radiationless  transition  1 Il>  •*  |II1>  or  Kp^  for  photochemical 
reaction.  and  f22^  correspond  to  the  absorption  system  (1I> 

•* — ►  1II>)  and  and  f^j^  correspond  to  the  emitting  system  (|III> 
-*  1I>.  Other  parameters  are  the  transition  frequencies  “21'  “31 
and  the  frequencies  u  and  u'  of  the  incident  and  eicitted 
radiation  respectively.  We  assvune  throughout  that 
=  Onr  aim  is  to  find  the  best  conditions  for  the 

enhancement  of  fluorescence  or  the  photochemical  decomposition 
rates . 


In  most  of  the  calculations  presented  here  we  use  the 

following  values  for  the  parameters:  a  =  200  A,  fjj^  =  f3j^  ®  0.1, 

=  3.48  eV,  «  =  3.46  eV  «  u ' ,  T-  =  r_  =  10®  sec 

K  =  10^  sec  and  the  incident  laser  power  of  lO^W/m^ .  We  call 
this  the  basic  parameter  set  and  mention  explicitly  only  the 
parameters  whose  value  differs  from  that  given  in  the  basic  set. 

For  a  small  sphere  the  reflection  tensor  R(u))  is 

R(u)  =  2a®(a  +  d)"®(e(u)  -  l)(e(w)  +  2)"^  (III.l) 

Since  we  assume  a  transition  dipole  perpendicular  to  the  surface 
and  along  the  2  axis  only  one  component  of  the  tensor  is  needed. 

In  Fig.  (2a)  we  plot  |R(u)|^  and  ll+R(u)l^  as  a  function  of 

(u,  for  several  values  of  d.  The  curves  show  that  the  local  field 

2 

enhancement  factor  |l+R(u)|  peaks  at  frequencies  which  are 


20 


slightly  lower  than  the  resonance  frequency  of  the  sphere  where 
I R{w) 1 ^  peaks) . 


For  a  small  sphere  .the  image  tensor  is 


G(u)  =I  {n+l)^a^”'*‘^(a+d)“^^"'^^^  (e{u)-l)  [e(u)  +  (n+l)/n]~^. 


n*l 


{III. 2) 


Since  we  assume  that  the  transition  dipole  is  perpendicular  to 
the  surface  and  lies  along  the  z  axis,  no  other  component  of  the 
tensor  is  needed.  While  the  reflection  tensor  R(u)  has  a 
resonance  at  a  frequency  given  by  Re  e(a))  =  -2.  (For  Ag  this  is 
at  u  =  3.49  eV,  if  the  dielectric  constants  of  Ref.  32  are  used), 
G(u)  has  an  infinite  niimber  of  resonances  whose  frequencies  are 
given  by  Re  6{u)  =  {n+l)/n.  Since  the  resonant  terms  are 
weighted  (in  Eq.  (III. 2))  by  the  factors  ( n*l )  ^  [a/ ( a+d)  ]  , 

their  relative  importance  depends  on  the  molecule-sphere  distance 
d.  If  d  is  very  small  only  the  resonance  corresponding  to  very 

7 

large  n  can  enhance  G((ij)  ;  this  gives  for  the  resonance  frequency 

Re  e(u)  =  -1  which  is  the  condition  for  a  flat  surface. 

Physically  this  makes  sense  since  if  d/a  <<  1  the  surface  appears 

flat  to  the  emitting  dipole.  At  large  values  of  d  the  dipole 

field  becomes  spatially  smooth  and  only  the  dipolar  resonance  n  = 

1  (i.e  Re  e  =  -2)  is  coupled  to  the  dipole.  Thus  we  expect  that 

the  line  width  is  resonantly  enhanced  through  ImG(u)  at 

frequencies  depending  on  d.  This  is  illustrated  in  Fig.  2b  where 

2 

we  plot  the  surface  induced  line  width  =  (e  f  ^ j^/2mu2,  )  ( P . 

P22)I™G(u)  as  a  function  of  u,  for  several  surface-molecule 
distances.  We  calculated  the  steady  state  values  of  -  P22' 
instead  of  using,  the  customary  assumption  that  p^^^  -  P22  -  1,  to 
test  whether  the  enhancement  of  the  local  field  causes  any  build 
up  in  022^  find  that  this  is  not  the  case.  The  calculation 
was  done  for  the  three  level  system,  that  is  for  K  ;«!  0 . 
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The  enhancement  factor  for  the  absorption  by  the  transition 
1 1>  |II>»  when  the  transition  |II>  -*■  |III>  is  turned  on  with 

the  rate  constant  K,  is  the  ratio  of  the  power  absorbed  in  the 
presence  of  the  sphere  P  ,  -to  the  power  absorbed  when  the  sphere 

Q 

is  absent  and  it  is  given  by 


R  1.  =  P  /  P 

abs  abs  abs  = 


(Pll“  P22^ 


ll+R{u)  +  ((r2+K)/2)^ 
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(III, 3) 


The  populations  and  are  given  by  the  steady  state  equa¬ 
tions  (11.33),  while  real  part  of  '>'21^'^^  given  by 

Eg.  (11.13).  The  enhcuicement  factor  for  f lu,orescence  is 

given  by  ,Sq.  (II. 36c).  In  Fig.  3  we  plot  and  R^,^  as 

functions  of  surface-molecule  distance.  Curve  (a)  has  r.  and  r_ 

9  —1  ^  ^ 

equal  to  10  sec  and  the  rate  coefficient  K  for  the  |I>  -►  |III> 

"1  -1 

transition  is  10  sec  .  Curve  (b)  has  an  increased  rate 
12  -1 

K  =  10  sec  ,  Since  this  increases  the  fluorescence  of  the 
molecule  at  the  surface  and  that  of  the  molecule  in  vacuum  we  do 
not  expect  a  dramatic  increase  of  the  enhancement  factor.  Some 
increase  is  achieved  because  the  rate  K  of  the  transition  |II>  -*■ 
|III>  competes  more  successfully  with  the  surface  quenching  of 
the  level  |II>  through  |II>  -*  1I>-  Curve  (c)  shows  the  effect  of 
increasing  which  Increases  the  fluorescence  rate  (the  scale 

for  curve  (c)  is  on  the  right  hand  side  of  Fig.  3)  and  the  rate 
of  absorption. 


The  enhancement  of  the  steady  state  population  (given 

in  Eq.  (11,33))  is  shown  in  Fig.  4,  for  the  same  parameters  as  in 

Fig.  3.  In  Fig.  5  we  show  the  relative  magnitude  of  the 

molecular  line  width  in  vacuum  with  respect  to  the  line  width 

in  the  presence  of  the  surface,  for  various  values  of  d  and  for  X 
12  - 1 

=  10  sec  .  Once  the  ratio  equals  1,  i.e.  once  r2 


exceeds 


substantially  the  surface  induced  width,  the  fluorescence  damping 
by  the  surface  is  no  longer  an  important  factor;  however,  the 
surface  still  can  enhance  fluorescence  through  enhanced 
absorption  {|I>  •*  |II>)  and  emission  -*•  |1>). 

In  Fig.  6  we  show  the  dependenc-s  of  the  absorption  and 
fluorescence  enhancement  factors  on  the  surface-molecule  distance 
for  several  choices  of  the  absorption  and  emission  frequencies  u 
and  w',  respectively.  For  all  cases,  K  =  10  sec  and  w  = 
and  u'  =  “3  2^*  i.e.  the  molecule  is  on  resonance  with  the  incident 
light  and  the  resonance  emission  is  detected.  Curve  (a)  has  u 
and  «'  close  to  each  other  and  in  resonance  with  the  sphere.  In 
curve  (b)  the  emission  frequency  is  off  resonance.  This 
diminishes  the  enhanced  emission  (i.e.  |l+R{(i)')|  )  but  it  also 
diminishes  the  rate  of  quenching  of  level  I11I>  by  the  surface; 
the  overall  effect  is  an  enhancement  of  the  fluorescence  for  the 
smaller  values  of  d.  The  other  curves  show  the  decrease  of 
and  as  both  frequencies  are  off  resonance. 

In  Fig.  7  we  show  the  enhancement  factor  R  for  the 

pc 

photochemical  rate,  (obtained  from  Eq.  (11.39)),  for  various 

values  of  the  rate  in  vacuum  K  .  The  enhancement  of  slow  rates 

pc 

is  fairly  marginal  at  all  distances,  essentially  because  surface 

quenching  of  state  |II>  is  more  efficient  than  its  depletion  by 

the  photochemical  process.  For  the  curve  (c),  having  the  highest 

value  of  quenching  process  prevails  at  low  distances  d, 

but  is  overcome  at  roughly  20  A.  The  reason  for  this  is  that  the 

2 

enhancement  factor  |1+R(«)I  decays  with  the  distance  slower  than 
the  quenching  rate;  therefore  at  some  intermediate  distance  it 
prevails  and  yields  an  enhancement  factor  of  50. 

In  Fig.  8  we  show  the  dependence  of  the  photochemical 

enhancement  factor  on  t.ne  Incident  power,  for  K  =  lO"^  sec”'^ 

3  8  10  2 

and  I,  at  10  ,  10  ,  and  10  W/m  .  Increased  power  diminishes 

me 

the  enhancement  factor,  essentially  because  in  vacuum  the  process 
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is  proportional  to  the  power.  The  slight  decrease  of  the 
enhancement  is  caused  by  the  fact  that  the  stimulated  emission 
from  level  |II>  is  enhanced  by  the  presence  of  the  sphere  and 
this  increases  the  radiative  energy  loss  from  level  |II>.  So 
while  the  absolute  photochemical  rate  goes  up  with  power,  the 
enhancement  by  the  sphere  goes  down.  For  high  powers  one  must 
also  be  concerned  with  the  effects  of  heating  and  multiphoton 
processes . 

In  concluding  we  emphasize  that  we  have  considered  here 

only  the  quenching  caused  by  energy  transfer,  as  given  by  local 

30 

electrodynamics.  The  non-local  effects  are  short  ranged  and 

are  important  only  at  the  distances  at  which  the  enhancement 

factors  are  very  low.  Therefore,  unless  one  is  specifically 

interested  in  monolayers  or  biiayers  the  non-local  effects  can  be 

disregarded.  Another  short  range  quenching  mechanism  is  the 

33 

charge  transfer  from  the  ejtcited  molecule  to  the  metal.  This 

33 

seems  to  be  the  predominant  mechanism  whenever  the  ionization 

potential  of  the  excited  state  exceeds  the  work  function  of  the 

34 

surface  by  less  than  approximately  1  eV.  We  have  made  no 
attempt  to  include  this  effect  here. 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  diagram  of  a  three  level  system,  u  and  u' 
are  the  incident  and  emitted  light  frequencies, 
respectively.  E^,  Ej*  and  E^  are  the  respective 
level  energies.  Solid  vertical  lines  indicate 
stimulated  transitions;  dashed  vertical  lines 
indicate  spontaneous  (radiative)  transitions;  the 
wavy  line  indicates  a  radiationless  transition. 

Figure  2(a).  |R(«)|^  (-.-•-)  and  |l+R(u)l^  { - — )  as  a  function 

of  frequency  u  for  various  d.  The  basic  parameter 
set  is  used  in  computation.  Curves  a,  b,  and  c  are 
for  d  -  10  A,  50  A,  and  150  A  respectively. 

Figure  2(b).  Surface  induced  width  as  a  function  of  frequency 

for  various  d.  The  basic  parameter  set  is  used  in 
computation.  Curves  a,  b,  and  c  refer  to  the  same 
situation  as  in  Figure  2(a). 

Figure  3.  Absorption  ( - )  and  fluorescence  ( - )  enhancement 

ratios,  Rgjjg  and  Rpj*  as  functions  of  distance  d 

from  the  surface.  We  use  the  basic  parameter  set, 

and:  (a)  K  *  10^^  sec~^;  (b)  «'  =  u.,  =  3.45  eV;  and 
10  -1 

(c)  =  10  sec  and  w'  as  in  curve  (b) .  The 

scale  for  curve  c  is  on  the  right  hand  side  of  the 
graph. 

Figure  4.  The  ratio  of  the  population  of  level  2  in  the 

presence  of  the  sphere  to  that  in  the  absence  of  the 
sphere,  as  a  function  of  d.  We  use  the  parameters  of 
curves  (a),  (b) ,  and  (c)  of  Fig.  3. 

Figure  5 .  The  ratio  of  the  natural  width  to  the  surface  induced 
width  as  a  function  of  the  natural  width.  Here  the 
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Figure  6. 


Figure  7 , 


Figure  8 , 


basic  parameter  set,  with  varying  is  used, 
a,  b,  and  c  are  for  d  »  20  A,  50  A,  and  100  A. 
respectively. 


Curves 


Absorption  { - )  and  fluorescence  (- 


— -)  enhancement 
ratios  as  ftinctions  of  the  distance  d.  The  various 
curves  are  obtained  by  using  the  basic  parameter  set 
with  the  following  changes: 
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The  enhancement  factor  (R__)  for  the  photochemical 

pc 

rate  as  a  function  of  distance  d.  The  basic 
parameter  set  with  the  replacement  of  K  by  K  is 
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